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Abstract 

We study predictions from perturbative Quantum Chromodynamics (QCD) for 
the process pp — >H — >77+X. In particular, we compare fully differential calculations 
at next-to-leading (NLO) and next-to-next-to-leading order (NNLO) in the strong 
coupling constant to the results obtained with the MC@NLO Monte Carlo (MC) 
generator, which combines QCD matrix elements at NLO with a parton shower al- 
gorithm. Estimates for the systematic uncertainties in the various predictions due to 
the choice of the renormalization scale and the parton distribution functions are given 
for the inclusive and accepted cross sections and for the corresponding acceptance 
corrections, obtained after applying standard selection and acceptance cuts. Further- 
more, we compare the distributions for the Higgs signal to those for the irreducible 
two-photon background, obtained with a NLO MC simulation. 
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1 Introduction 



One of the main goals of the future Large Hadron Collider (LHC) at CERN is the search 
for the Higgs boson. From theoretical considerations [1] it is concluded that the Higgs 
mass should be below about one TeV/c 2 . In addition, precision measurements of electro- 
weak observables [2| indicate that the low- mass region (approximately below 200 GeV/c 2 ) 
is preferred. If the Higgs boson mass is smaller than twice the W mass, its decay mode 
into a pair of photons will be among the most important and distinct channels for the 
Higgs search. Since this decay can only proceed via loop corrections, the corresponding 
branching ratio is less than one percent. However, it provides the experimentally cleanest 
signature jHllHE]- 

Obviously, one important aspect of the preparations for the upcoming LHC Higgs 
searches is to study the predictions for the distributions of kinematic variables for the 
signal, pp — > H — > 77 + X, and the background processes (eg. pp — > 77 + X), obtained 
from perturbative Quantum Chromo dynamics (QCD). First, it is known that the inclusive 
Higgs cross section gets large perturbative corrections at next-to-leading order (NLO) 
001313 and next-to-next-to-leading order (NNLO) flU EJ E2J Q31 EE] , so it is likely that 
the dynamics of the process is different at various orders in perturbation theory. That is, the 
QCD corrections to the production process change the kinematics of the Higgs boson and 
lead to modifications in the kinematics of the final-state photons. It is worth noting here 
that very recently the leading threshold-enhanced third-order (N 3 LO) corrections to the 
total Higgs cross sections have been derived showing that the perturbative expansion 
stabilizes. 

Second, the detailed knowledge of differential cross sections will help to optimize the 
search strategy and thus improve the signal over background ratio. In the past, such 
studies have been performed using leading order (LO) Monte Carlo generators, which 
usually include a parton shower and a hadronization model and are finally combined with 
a full detector simulation. Higher order QCD predictions in general were only available 
for the inclusive cross sections or for distributions of some kinematic variable of the Higgs 
boson, for example its transverse momentum, [TBI IT7] . Therefore, the X-factors applied 
to the LO predictions to correct for the missing higher order terms would not take account 
of the analysis selection and acceptance cuts, which might reduce the allowed phase space 
in a non-trivial way. An attempt to improve on this consists in a re-weighting procedure 
where the LO plus parton shower results are re-scaled in order to reproduce the higher 
order prediction for a particular kinematic variable, such as the Higgs pp- Taking account 
of these weights and of the efficiencies obtained from the MC simulation, an effective 
.fT-factor can be computed which incorporates the selection and acceptance factors [18]. 

Recently, new calculations for the process pp — >H —^77 + X beyond NLO have be- 
come available, which now allow for an analysis of acceptance corrections and the related 
systematic uncertainties at higher order in perturbative QCD. In particular, the fully dif- 
ferential cross section at NNLO has been computed for the first time [T§1 I20j. giving the 
exact expressions (up to this order of perturbative QCD) for Higgs production and its 
two-photon decay with up to two additional massless partons (quark or gluons). Therefore 
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this new semi-analytical calculation allows to apply realistic selection and isolation cuts to 
the Higgs decay products and the additional partons in the final state, as well as to study 
the distribution of relevant kinematic variables, such as the rapidity and pt of the Higgs 
and the photons, at NNLO accuracy. 

An alternative approach consists in the combination of QCD matrix elements at NLO 
and a parton shower algorithm, therefore taking into account also large logarithms due 
to soft and collinear radiation at all orders in the strong coupling constant. This is im- 
plemented in the Monte Carlo (MC) generator MC@NLO [2H 122], combined with the 
HERWIG MC program [231211, which in addition handles the hadronization process and 
particle decays. Several aspects render such an approach interesting for phenomenological 
studies and experimental analyses. In contrast to a pure parton level calculation, such 
a MC generator allows to investigate distributions at the particle (hadrons and leptons) 
level, which is closer to the actually measured final state than partons. For example, the 
detailed internal structure of jets is modeled rather well, studies of lepton isolation from 
jets are expected to be more realistic and it is possible to simulate also the effects of the 
underlying event and multiple interactions. Finally, it is rather straightforward to feed the 
output of this generator into a full detector simulation, which is necessary for the experi- 
menter in order to develop optimised data analyses which rely on a good approximation of 
the underlying QCD dynamics. There is already a considerable list of processes available 
in MC@NLO, among these pp — >H — >77 + X, and further additions can be expected in 
the future. 

This paper provides a first detailed comparison of these newly available predictions. 
Taking the various calculations, we apply standard cuts to relevant kinematic observables 
in order to reproduce in a realistic manner the basic aspects of the Higgs search in the two- 
photon channel at LHC. As a result, we are able to compare the predictions for the inclusive 
and accepted cross sections, as well as for the corresponding acceptance corrections. These 
comparisons are done as a function of the Higgs mass and the renormalization/factorization 
scales. In addition, we study the dependence of the MC@NLO predictions on the choice 
of the parton distribution functions (pdf). Regarding the latter, it can be interesting to 
compare the predictions when combining MC@NLO with pdfs extracted at NLO as well as 
at NNLO accuracy. Although via the combination of NLO matrix elements with the initial 
state parton shower effectively also terms beyond NLO are generated for the perturbative 
part of the scattering process, the use of NNLO pdfs in MC@NLO should be regarded as 
equivalent to that of NLO pdfs, since, formally, the accuracy of the results is NLO in both 
cases. 

Interesting kinematic observables of the two-photon final state, as proposed in |25 j and 
[20], are computed for the signal as well as for the irreducible two-photon background, 
pp ^77 + X. The latter is known at NLO accuracy [23 120] ■ It is shown that variables 
such as the average photon transverse momentum or the rapidity difference of the two 
photons could contribute to the discrimination of the signal from the background. 

It is worth noting that the entire study is done at the phenomenological level, i.e. no 
attempt is made here to simulate any detector effects beyond the basic selection criteria, 
such as acceptance cuts in photon rapidity and transverse momentum, as well as isolation 
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requirements. We thus emphasize that this paper does not provide a re-evaluation of 
the Higgs discovery potential at LHC, but rather intends to show the level of theoretical 
accuracy, which is achievable with the tools that have become available recently, as well 
as to indicate where more detailed studies, including full detector simulation and the most 
up-to-date perturbative QCD predictions, might be worthwhile pursuing. 

The paper is organized as follows: In Section |2] we list the details on the computer 
programs and the relevant parameter settings which were used to perform this study. In 
Section 01 the basic event selection cuts are described. The resulting cross sections before 
and after this selection are discussed in Section H] and the distributions of some kinematic 
variables are presented in Sectional The comparison to the main background is given in 
Section |H1 Finally, Section [7| contains a summary of the results. 

2 Computer programs and parameter settings 

In order to carry out this study we have used several different software packages, which 
implement the various signal and background predictions up to a certain level of approxi- 
mation in perturbation theory. The employed programs are : 

• FEHIP [20] (Version 1.0) 

This code is used to calculate the fully differential cross section for Higgs boson 
production through gluon fusion at NLO and NNLO QCD, in the narrow-width 
approximation. The gluon fusion channel via a top quark loop contributes a fraction 
of about 80% to the total Higgs production cross section in the relevant low-mass 
region. The calculation has been carried out in the limit of an infinitely heavy top 
quark, in which the Higgs boson coupling to gluons is point-like. However, the results 
are re-scaled using the ratio of the LO cross section predictions, obtained with the 
exact top mass dependence and in the heavy top mass limit. As discussed in (20], 
this gives an excellent approximation. 

For the consistent convolution of the partonic cross sections at NLO (NNLO) with 
parton distribution functions, the MRST01 NLO (NNLO) j2Z] pdf sets are used. The 
so-called mode 1 is employed as the default mode for pdf evolution. For the moment 
being, no other pdf set is available with this program and we have made no attempt 
to include further sets. 

The code returns the four-momenta of the Higgs boson and the final state particles, 
namely two photons and up to two additional massless partons, which can be used 
to evaluate an analysis function for every event. This analysis function is then con- 
voluted with the squared matrix elements of the corresponding final state in order to 
finally obtain a cross section. Many of the numerical results, in particular at NNLO, 
were obtained from |28*1 120j. Concerning statistical uncertainties, the NNLO results 
for the cross sections (acceptance corrections) are accurate to 1% (2%), while the 
NLO numbers have been computed up to an accuracy of 0.1%. 
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In order to speed up the computation, for the Monte Carlo integration step a parallel 
VEGAS algorithm has been used [23 EES- The available C++ code has been slightly 
modified in order to meet the specific needs of FEHIP. The code has been parallelized 
using the Standard Message Passing Interface MPI [3*T] . 

• MC@NLO [HI 122] (Version 2.31) 

As described in the introduction, this MC generator consistently incorporates NLO 
QCD matrix elements into a parton shower framework. In the purely inclusive case 
the code returns the NLO cross section, whereas the effects of the parton shower 
become only apparent when cuts are applied to the final state. In order to eliminate 
a problem with this version of the code when using renormalization and factorization 
scales that differ from the default choice {fi = mn), a correction has been applied 

MC@NLO is used together with the HERWIG event generator [2Hll2Ij, which handles 
the simulation of the parton shower, the hadronization process and possible particle 
decays. The exact top mass dependence is retained at the Born level. We have 
generated events using three different pdf sets. In addition to the two MRST01 sets 
mentioned above, we have also used the CTEQ6M set, which has been extracted 
at NLO accuracy. 

• DIPHOX (Version 1.2) 

The background process pp — ^77 + A 7 " is computed with this program, which provides 
a NLO QCD prediction at parton level only. The (unphysical) phase-space slicing 
parameter pr m j2Hj was set to 0.05 GeV/c. An additional contribution to this process, 
which is not included in DIPHOX, namely the two-loop correction to the gluon fusion 
subprocess, has been calculated in There it is shown that this contribution is 
modest, around 10% or smaller for the mass range of interest in the low-mass Higgs 
search. 

• HDECAY |3U (Version 3.101) 

This is used to calculate the branching ratio of the Higgs decay to two photons, 
resulting in BR(H ^77) = 0.22 %. 

• HIGLU |3S] (Version 2.102) 

The program HIGLU allows to calculate the inclusive cross section for Higgs boson 
production through gluon fusion at NLO, with the exact top mass dependence at 
both LO and NLO. 

Throughout this analysis a top mass of m t0 p = 175 GeV/c 2 has been assumed. When- 
ever we refer to 'the scale' /1, it means /i = /ir = /if, i.e. the factorization and renormal- 
ization scales are set to the same value. 
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3 Event selection 



In this study we attempt to simulate in an appropriate manner the planned Higgs searches 
at LHC in the two-photon channel, using only phenomenological, i.e. parton and/or parti- 
cle level, predictions without any detailed detector simulation. This is achieved by taking 
into account the geometrical limitations of a typical LHC detector and the foreseen trig- 
ger requirements. Furthermore, in order to reduce the large irreducible background the 
following cuts are applied to the final state photons : 

• The absolute value of the pseudorapidity |% 12 | of each of the two photons must be 
smaller than 2.5. This corresponds to the coverage in polar angle of a detector such 
as ATLAS [H] or CMS [37]. 

• One of the two photons must have a transverse energy above 40 GeV, while the other 
photon is required to have a transverse energy above 25 GeV. 

• In order to exclude events with considerable hadronic activity around the final state 
photons, as for example can be expected in background events with quark-photon 
fragmentation, the photons must be isolated, i.e. the additional energy in a cone 
AR = y/ Ai] 2 + Aip 2 < 0.4 around each photon must not exceed 15 GeV. 

This selection is basically the same as the 'standard cone criterion' used in [23]. For 
reasons given in that reference (e.g. calorimeter granularity and finite photon shower size), 
we did not investigate the smooth cone isolation criterion jHEj- In the following, we will 
refer to these selection cuts as the 'standard cuts' and call the accepted photon with the 
highest (second highest) transverse energy photon 1 (photon 2). 

4 Inclusive and accepted cross sections 

In the following section we study the predictions at various orders in perturbative QCD 
for the inclusive pp — >H — >77 + X cross section a inC) as well as for the accepted cross 
section cr acc found after the standard cuts on the final state photons have been applied. 
The comparison is made for several different choices of parameters, such as the Higgs mass, 
the scale or the pdf set. For the moment we simply assume that the branching ratio of the 
Higgs decay into two photons is one. 

In a first step, we fix the Higgs mass to m H = 120 GeV/c 2 , the factorization scale /ip 
and renormalization scale /zr are set to mn/2 and the pdfs are taken from the MRST01 
set, using the NLO (NNLO) pdfs in combination with MC@NLO and FEHIP at NLO 
(NNLO). As discussed in detail in Ref. [20] , the choice /if = /iR = mu/2 is motivated 
by the observation of an improved convergence of the perturbation series and the very 
good agreement with the threshold-resummed results for the Higgs hadroproduction cross 
section |55] . This is further confirmed by the recent calculation of the threshold-enhanced 
third-order corrections [To] . 
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The results can be found in Tabled Whereas with MC@NLO we give the acceptances 
for each single cut as well as for groups of cuts, in the case of FEHIP we had to restrict 
ourselves to the acceptance after all cuts, since this NNLO code is very cpu-intensive. 
Regarding the inclusive cross section, a- mc , we observe that the MC@NLO prediction falls 
short of the NLO prediction of FEHIP by about 3.3%. By construction MC@NLO should 
return the exact NLO cross section if no cuts are applied. Thus the observed difference has 
to be explained by the different treatment of the top mass dependence. FEHIP uses the 
infinite top mass limit in all orders of the calculation and then re-scales the result using 
the ratio of the exact and approximated (mt op — > oo) LO cross sections, while MC@NLO 
keeps the exact top mass dependence at LO, but uses the m top — > oo approximation at 
NLO. For comparison, the result obtained with HIGLU which retains the exact top mass 
dependence at LO and NLO is a inc = 44.39 ± 0.01 pb. Thus FEHIP returns a 0.9% higher 
and MC@NLO a 2.5% lower NLO cross section with respect to HIGLU. 

The differences of the NLO results with respect to the NNLO result amount to 13% 
(MC@NLO) and 9% (FEHIP, NLO). After the standard selection, these differences are re- 
duced to about 10% and 6%, respectively. The acceptance corrections for m H = 120 GeV/ c 2 , 
defined as the ratio between the accepted and the inclusive cross sections, cr acc /cr inc , agree 
to within 2% (absolute). 

It turns out that the NNLO calculation predicts a slightly stronger reduction in the 
accepted cross section than the other two approaches. With the MC@NLO generator 
we have also studied the effect of the individual cuts. We find that the reduction in cross 
section is mainly due to the angular and momentum cuts, whereas the isolation requirement 
has only a very minor impact. 

4.1 Choice of pdfs in MC@NLO 

As a first variation with respect to the parameter choices given above, we investigate 
the dependence on the pdfs used together with MC@NLO. Since MC@NLO includes the 
leading effects of soft and collinear radiation to all orders, thus beyond NLO, it is interesting 
to combine it with the NNLO pdf set of the MRST01 fits, which has been extracted using 
a combination of NLO and NNLO cross sections. Finally, we also take CTEQ6M (NLO) 
as input, in order to study an alternative set of pdfs, obtained under different assumptions, 
from a somewhat different data set and thus with different theoretical and experimental 
systematic uncertainties. 

The results are listed in Tabled Compared to MRST01 NLO, the results found with 
the MRST01 NNLO set, for both cr inc and cr acc , differ quite considerably, up to 16%, whereas 
the cross sections obtained with the NLO sets from the MRST01 and CTEQ6M fits agree 
to better than 4%. The considerably lower cross sections obtained with the NNLO set are 
interpreted as arising from the smaller gluon density and strong coupling used in this set. 
It is worth noting that in the NNLO calculation, as implemented in FEHIP, the NNLO 
corrections in the matrix elements, which also include additional scale dependent terms, 
counter the effect of the NNLO evolution of the pdfs. 

Interestingly, when looking at the acceptance, we find agreement for all three choices 
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Generator 


MC@NLO 


FEHIP NLO 


FEHIP NNLO 


PDF set 


MRST01 NLO 


MRST01 NNLO 


CTEQ6M 


MRST01 NLO 


MRST01 NNLO 


0"inc [Pb] 


43.30 


37.67 


44.32 


44.77 


48.94 





Acceptances (a acc /a inc ) 


of single groups o 


f cuts 


PT-CUtS 


80.2 % 


80.3 % 




80.1 % 






r^-cuts 


83.2 % 


83.0 % 




83.8 % 






isolation 


99.7 % 


99.7 % 




99.8 % 







Acceptances (o" acc /cri nc ) of two groups of cuts 


Pt- and r?-cuts 


63.6 % 


63.2 % 


64.3 % 






PT-cuts and isolation 


79.6 % 


79.7 % 


79.6 % 






??-cuts and isolation 


81.5 % 


81.1 % 


82.1 % 







All three groups of cuts 


Acceptance 


63.0 % 


62.6 % 


63.7 % 


63.4 % 


61.4 % 


Oacc [pb] 


27.30 


23.60 


28.25 


28.35 


30.07 



Table 1: Cross sections and acceptance corrections for a Higgs mass of m H = 120 GeV/c 2 , obtained with MC@NLO and 
FEHIP (NLO and NNLO) for different cuts and pdf sets. The renormalization and factorization scale are set to ji — m H /2. 
The branching ratio H —^77 is assumed to be one. The statistical uncertainties are at the one (two) per-cent level for the 
NNLO cross sections (acceptance) and negligible for the other cases. 



within 1.1%, absolute; when also considering the FEHIP results the span between the 
largest and smallest acceptance increases to 2.3%. Half of this span could be defined as a 
1.8% 'uncertainty' relative to the average of the largest and smallest values. Calculating the 
same relative uncertainty for the accepted cross section, we find 4.7%. Thus the acceptance 
correction is less sensitive with respect to the choice of pdfs and the approximations made 
by the perturbative calculation than the absolute cross section. 

4.2 Variation of the Higgs mass 

In TableElthe results are given for the inclusive and the accepted cross sections for different 
Higgs masses, obtained with MC@NLO and FEHIP (NNLO). The renormalization and 
factorization scales are chosen to be fi = mn/2. 

In the upper part of the table the inclusive cross sections a- inc for MC@NLO (MRST01 
NLO) and the NNLO results obtained with FEHIP (MRST01 NNLO) are compared. It 
can be observed that the inclusive cross sections obtained with MC@NLO are 16% smaller 
than the NNLO results for the full Higgs mass range investigated here. In the second part 
of TableElthe MC@NLO predictions for the acceptance corrections, after the various cuts, 
are given. As stated above, the isolation cut has the smallest effect on the acceptance 
corrections. The increased acceptance of the cuts on pr arid r\ when increasing the Higgs 
mass arises from the higher transverse momenta expected in this region. Therefore also the 
total acceptance increases with larger mn- The results after cuts, obtained with MC@NLO 
and FEHIP, agree to within 11%. This reduced difference compared to the inclusive case 
is due to the slightly larger acceptance correction found with MC@NLO. Remarkably, the 
acceptance corrections for MC@NLO and FEHIP (NNLO) agree within 3.2% (absolute) 
over the whole investigated mass range. 

The results are shown in Fig. where the inclusive and accepted cross section, as well 
as the acceptance correction, are plotted as a function of mg. It is worth noting that the 
statistical uncertainty of the NNLO results for the acceptance correction is at the 2% level. 

4.3 Variation of the renormalization and factorization scale fi 

Both the renormalization and factorization scales are un-physical variables which are in- 
troduced in the regularization step of the renormalization and/or factorization procedure, 
which is needed to deal with the divergent integrals encountered in the perturbative cal- 
culations. An approximation for a perturbatively calculated observable, obtained to some 
order in the coupling constant, does not depend on these scales up to the calculated order, 
but a dependence appears at the next-higher order. It is thus interesting to investigate 
the remaining dependence of a given higher-order calculation. The range over which these 
scales are varied, and the usually quoted systematic uncertainties corresponding to the 
spread in the results when varying the scales, are somehow arbitrary and subject to con- 
troversy. For this analysis we have studied scale variations between mn/4 and 2mn, in 
both MC@NLO and FEHIP. It is worth noting that other settings of the scales are possi- 
ble in MC@NLO, for example, including some pi-dependence. Our approach is convenient 
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Generator 


MC@NLO 


NNLO 


MC@NLO 


NNLO 


MCONLO 


NNLO 


MC@NLO 


NNLO 


MC@NLO 


NNLO 


m H [GeV] 


110 


120 


130 


140 


150 




50.65 


56.82 


43.30 48.94 


37.45 


42.54 


32.66 


37.37 


28.70 


33.10 




Acceptances (o" acc /<7 inc ) of single groups of cuts 


PT-cuts 

?7-cuts 

isolation 


75.9 % 
82.4 % 
99.7 % 




80.2 % 
83.2 % 
99.7 % 




83.3 % 
84.0 % 
99.7 % 




85.9 % 
84.8 % 
99.7 % 




87.8 % 
85.3 % 
99.7 % 





Acceptances (cr acc /o~i nc ) of two groups of cuts 


Pj- and //-cuts 
PT-cuts and isolation 
//-cuts and isolation 


59.9 % 
75.4 % 
80.7 % 




63.6 % 
79.6 % 
81.5 % 




66.3 % 
82.7 % 
82.1 % 




68.6 % 
85.2 % 
82.8 % 




70.5 % 
87.1 % 
83.3 % 





All three groups of cuts 


Acceptance 


59.5 % 


57.3 % 


63.0 % 


61.4 % 


65.7 % 


63.7 % 


67.9 % 


64.7 % 


69.7 % 


67.4 % 


CTacc [Pb] 


30.11 


32.56 


27.30 


30.07 


24.60 


27.08 


22.19 


24.17 


20.01 


22.30 



Table 2: Acceptance corrections for the different cuts and the different Higgs masses, as obtained with MC@NLO and 
FEHIP (NNLO). The used pdfs are MRST01 NLO for MC@NLO and MRST01 NNLO for FEHIP. The renormalization 
and factorization scale are set to /j, — toh/2. The statistical uncertainties are at the one (two) per-cent level for the NNLO 
cross sections (acceptance) and negligible for the other cases. 
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Figure 1: Inclusive and accepted cross sections as a function of the Higgs mass, mn, 
obtained with MC@NLO and FEHIP (at NLO and NNLO). 

for a direct comparison of MC@NLO and FEHIP. 

The NNLO results when choosing either /x = mn/2 or /i = 2mn are shown in Fig. El 
As can be observed in Fig. El the NLO cross sections start to increase very rapidly when 
going to very small scales, whereas the scale dependence at NNLO is considerably flatter. 
The difference between the NLO predictions from FEHIP and MC@NLO can be explained 
by the different treatment of the top mass dependence. 

For a Higgs mass of m H = 120 GeV/c 2 , as used in Fig. El it appears that a range of 
scales between mn/2 and 2mn might be reasonable for evaluating a possible systematic 
uncertainty because of this effect. Table El summarizes the MC@NLO results for the 
different scale choices. Despite the rather large variation of the cross sections when varying 
/j,, the acceptance correction is remarkably stable to within 0.5% (absolute), even when 
going to scales as low as mn/4. A similar stability of the acceptance correction has been 
found for single W production at LHC [40J. 
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Figure 2: Inclusive and accepted cross sections at NNLO as a function of the Higgs mass, 
?71h, for two different scale choices. 

5 Differential cross sections 

After having calculated the acceptance for a rather standard set of basic selection cuts, it 
is interesting to investigate further kinematic observables which might help to discriminate 
the signal from backgrounds. In order to be able to optimize such an additional selection, it 
is necessary to have a good understanding of the relevant kinematic distributions. Here we 
concentrate on two observables, which have also been discussed in Refs. |2*B^ |2*U] , namely 
the distributions of the mean of the transverse momenta of the two final state photons 
Pt = (Pt + Pt ) an d the difference of their pseudo-rapidities Y* = |^ 71 — rf l2 \ /2. As 
shown in [23], the latter is interesting since a similar distribution for the prompt photon 
background is flatter. In the following we focus on a Higgs mass of mn = 120 GeV/c 2 
and set the renormalization and factorization scale to fi = mn/2. All the distributions are 
obtained after applying the selection cuts as described in Sectional 

Figure |U shows the p™ (left) and Y* (right) distributions obtained with MC@NLO when 
using two different pdf sets, MRST01 NLO and CTEQ6M NLO. We observe that in the 
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Figure 3: The accepted cross section cr acc predicted by MC@NLO and FEHIP (NLO and 
NNLO) as a function of the renormalization and factorization scale, /i, for a Higgs mass of 
m H = 120 GeV/c 2 . 



interesting region 30 GeV/c < < 100 GeV/c the differential cross section found with 
CTEQ6M exceeds that obtained with the MRST01 NLO set by up to 5%, in accordance 
with the differences seen for the inclusive and accepted cross sections listed in Table [TJ 
Obviously, the same feature can be found in the Y* distribution. However, in this case it 
is remarkable that the ratio of the two predictions is flat over the whole Y* range, which 
indicates that the shape of this distribution is rather insensitive to the choice of the pdf 
set. 

The same distributions are reproduced in Fig. El where now the CTEQ6M pdfs have 
been replaced by those from the MRST01 NNLO set. Similarly to the observations above, 
the shape of the two histograms differ only slightly, in particular in the case of Y*, with 
the results based on the NNLO pdfs being lower than those from the NLO set. This has 
been discussed previously, cf. Section 14.11 and Table HJ Thus the assumption is confirmed 
that the choice of the pdf set does not strongly affect the shape of these distributions. 

Finally, in Fig. |U] we compare the distributions obtained with MC@NLO and with 
FEHIP at NNLO. Concerning the variable p™, we observe that in the peak region the shapes 
of the distributions differ rather strongly, with the maximum shifted to lower values when 
going from MC@NLO to NNLO. The appearance of these large perturbative corrections 
close to the kinematic boundary has also been discussed in (201 • From the upper tail, thus 
far beyond this boundary, it is evident that at larger photon momenta the cross section is 
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scale /i = /if = /iR 


H = m h /4 


H = mh/3 


/i = m h /2 


/i = mh 


/i = 2mh 


0"inc [pb] 


51.73 


47.85 


43.30 


37.16 


32.42 




Acceptances (cr acc /cr inc ) of single groups of cuts 


Px-cuts 

77-cuts 

isolation 


80.2 % 
83.2 % 
99.7 % 


80.2 % 
83.2 % 
99.7 % 


80.2 % 
83.2 % 
99.7 % 


80.2 % 
83.4 % 
99.8 % 


80.1 % 
83.5 % 
99.8 % 



Acceptances (o"acc/0"inc) of two groups of cuts 


Pt- and 77-cuts 
PT-cuts and isolation 
77-cuts and isolation 


63.5 % 

79.6 % 
81.3 % 


63.4 % 
79.6 % 
81.3 % 


63.6 % 
79.6 % 
81.5 % 


63.8 % 
79.7 % 
81.7 % 


63.8 % 
79.6 % 

81.9 % 




All three groups of cuts 


Acceptance 


62.9 % 


62.8 % 


63.0 % 


63.3 % 


63.4 % 


Oacc [Pb] 


32.52 


30.05 


27.30 


23.53 


20.54 



Table 3: Acceptance corrections for different cuts and scale choices, as obtained with 
MC@NLO. The Higgs mass is set to m H = 120 GeV/c 2 . 

higher at NNLO with respect to MC@NLO. Regarding the Y* distribution it is particularly 
interesting to observe that the step in the NNLO prediction around Y* = 0.9 — 1.0 does 
not appear for MC@NLO. Such a step is typical for fixed order calculations when a phase- 
space boundary varies at different orders in the perturbation series, however, it can be 
removed by the inclusion of soft-gluon resummation to all orders jH]. In the case of Y*, 
the leading order cross section is zero above Y* = 0.96. 

6 Comparison to the irreducible background 

After having analyzed the signal cross sections as functions of the measurable final state 
photon momenta, it is natural to compare these differential results to those for the irre- 
ducible background, namely prompt photon production, pp ^77 +X. The latter have been 
computed with DIPHOX|Z5| at NLO (cf. Section HJ), with MRST01 NLO as pdf set. The 
renormalization and factorization scales are fixed to \x = mn/2 for the signal, computed 
with MC@NLO, and \x = m 77 /2 for the background, where m 77 is the invariant mass of 
the two final-state photons. For the photon fragmentation we use the NLO fragmentation 
functions (set I) from |4*2*] . 

Figure [7| compares the normalized distributions of (left) and Y* (right) for signal 
and background. For the background, only events with 118 GeV/c 2 < m 77 < 122 GeV/c 2 
are included, corresponding to a mass-window of ±2 GeV/c 2 around the Higgs mass jtih = 
120 GeV/c 2 . As expected, the shapes differ considerably. In the case of p™, the mean is 
clearly shifted towards lower values for the background. The Y* distribution offers an even 
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Figure 4: The distributions of the kinematic variables (left) and Y* (right), obtained 
with MCONLO for the two pdf sets MRST01 NLO and CTEQ6M NLO. The standard 
cuts as indicated in the plots have been applied. 

better handle for the signal-background discrimination, since the background is rather flat, 
whereas the signal is peaked for small values of Y * . 

6.1 Statistical significance 

In Table0]we compare the signal cross section obtained with MC@NLO to the background 
cross section computed with DIPHOX. Here, <7i nc refers to the inclusive pp — *-H —^77 + X 
cross section before any cut. In addition to the acceptance correction found with MC@NLO 
(cf. Section HJ), we also account for a photon reconstruction efficiency of 57 % as was done 
in Ref. corresponding to a combination of 81% for 7/jet identification per photon and 
87% for fiducial cuts. The branching ratio for the Higgs decay channel into two photons 
is calculated with HDECAY. Altogether this gives an effective signal cross section of 
cx e ff = 0.034 pb. If we assume an integrated luminosity of 30 fb _1 , corresponding to about 
three years of LHC running at low luminosity, we find a total number of S = 1027 signal 
events. 

In order to calculate the number of background events, we only consider events with a 
total invariant mass of the two final-state photons between 118 and 122 GeV/c 2 . In this 
case, the cross section after cuts amounts to cr acc = 0.944 pb, as shown in the right column 
of TableEJ In addition to the efficiency factor, we include a 20 % reducible background [2o] . 
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Figure 5: The distributions of the kinematic variables (left) and Y* (right), obtained 
with MC@NLO for the two pdf sets MRST01 NLO and MRST01 NNLO. The standard 
cuts as indicated in the plots have been applied. 



This leads to an effective background cross section of a e g = 0.631 pb and to a total number 
of B = 18928 background events. Thus, we end up with a significance of S/\/B = 7.5. 

A few remarks are in place here. First, this is the significance obtained only after ap- 
plying the basic selection cuts given in Section |3J without any further signal-background 
discrimination as suggested, eg. by the observations from Fig. [7| Second, DIPHOX does 
not contain the two-loop correction to the gluon fusion subprocess of pp — >77 + X, which 
has been calculated in However, its inclusion would only lead to a 5% reduction of 
the statistical significance given in Table |3J Furthermore, it is worth noting that the back- 
ground has been calculated at NLO, while the signal obtained with MC@NLO includes 
corrections beyond that order. We think that a calculation of the prompt photon back- 
ground to such an approximation as can be obtained with the MC@NLO approach would 
be highly desirable, both from the phenomenological and experimental point of view. 



7 Summary 

We have analyzed perturbative QCD predictions at and beyond next-to-leading order accu- 
racy for the Higgs production and decay in the channel pp ^77 + X at LHC energies. 
In particular, for the first time it has become possible to compare the results from the 
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Figure 6: The distributions of the kinematic variables (left) and Y* (right), obtained 
with MC@NLO (pdf set MRST01 NLO) and FEHIP (NNLO). The standard cuts as indi- 
cated in the plots have been applied. 

MC@NLO Monte Carlo generator to a full NNLO prediction, not only for inclusive cross 
sections, but also for the accepted cross sections obtained after applying standard selection 
cuts. 

Whereas the inclusive and accepted cross sections at NLO and NNLO differ by about 
10%, it turns out that the acceptance corrections agree to within 2%. The latter are 
also very stable under variations of the renormalization and factorization scale and rather 
insensitive to the choice of the parton distribution functions. This is not the case when 
looking at the cross section themselves, in particular at NLO approximation. 

We have also analyzed distributions of kinematic observables constructed from the 
final-state photon momenta, namely the average of the two photon transverse momenta 
and their rapidity difference. Some differences in the shapes of the distributions have been 
found between MC@NLO and the exact NNLO calculation. More importantly, it has been 
shown that these distributions should allow for a good discrimination between the Higgs 
signal and the irreducible prompt photon background. 
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Figure 7: Normalized distributions of Y* and for signal and background. The standard 
cuts as indicated in the plots have been applied. 
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